Temperature dependence of spinon and holon excitations in one-dimensional Mott 

insulators 
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Motivated by the recent angle-resolved photoemission spectroscopy (ARPES) measurements on 
one-dimensional Mott insulators, SrCu02 and Nao.96V20s, we examine the single-particle spectral 
weight of the one- dimensional (ID) Hubbard model at half-filling. We are particularly interested 
in the temperature dependence of the spinon and holon excitations. For this reason, we have 
performed the dynamical density matrix renormalization group and determinantal quantum Monte 
Carlo (QMC) calculations for the single-particle spectral weight of the ID Hubbard model. In the 
QMC data, the spinon and holon branches become observable at temperatures where the short-range 
antiferromagnetic correlations develop. At these temperatures, the spinon branch grows rapidly. 
In the light of the numerical results, we discuss the spinon and holon branches observed by the 
ARPES experiments on SrCu02- These numerical results are also in agreement with the temperature 
dependence of the ARPES results on Nao.9eV205. 

PACS numbers: 71.10.Fd, 71.45.-d, 79.60.-i, 74.72.Jt 



I. INTRODUCTION 

One of the central issues in one-dimensional (ID) cor- 
related electron systems is the spin-charge separation Q . 
In these systems, the spin and charge degrees of free- 
dom of electrons are decoupled into their collective exci- 
tations, 'spinon' and 'holon' Since the single-particle 
excitations are not quasiparticles, it is expected that 
these excitations give rise to completely nontrivial phys- 
ical properties. 

Experimentally, angle-resolved photoemission spec- 
troscopy (ARPES) gives direct information on the single- 
particle excitation spectra. High-quality ARPES mea- 
surements under various physical conditions enable us to 
understand the collective excitations in ID compounds. 
Recently, such ARPES mearurements have been per- 
formed on a ID Mott insulator SrCu02, where the 
spinon and holon branches have been unambiguously ob- 
served 0- In this work, the spectral weight from the 
main valence band with 02p character is suppressed due 
to cross section effects, and thus the two branches become 
observable in the low binding-energy region. The band 
dispersions of the spinon and holon branches are in good 
agreement with those predicted in a spin-charge sepa- 
rated model 4J. Furthermore, careful lineshape analysis 
reveals that the peak height of the holon branch is smaller 
than that of the spinon branch, and the full widths at 
half-maximum of the spinon and holon branches are es- 
timated to be w 0.7 eV and w 0.5 eV, respectively. 

Temperature dependent ARPES studies are also im- 
portant, because finite-temperature effects on the single- 
particle excitation spectra of ID Mott insulators are not 
due to simple thermal broadening. Since the single- 
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particle excitation spectra are given by the convolution 
of the spinon and holon Green's functions, the finite- 
temperature effects reflect their collective nature, not the 
normal Fermi distribution of the quasiparticles. A tem- 
perature dependent ARPES study was carried out on 
Nao.96V205 0- In this work, the spectral weight re- 
distribution from higher to lower binding-energy region 
was observed when the temperature was decreased from 
300 K to 120 K. The redistribution occured on the scale 
of I eV, which was 100 times larger than the tempera- 
ture change. The ARPES data were consistent with the 
finite-temperature exact-diagonalization calculations for 
the ID t-J model at half- filling, and it was shown that the 
spin-charge separation picture is valid for Nao.96V2 0s. 

In the high-energy ARPES measurements on 
SrCu02 it is quite significant to detect the spinon 
and holon branches directly in the ARPES spectrum, 
because the ARPES data enable us to examine the 
spectral weights and the lifetimes of the spinon and 
holon in the compound. The temperature dependent 
ARPES measurement on Nao.geV^Os tells us how the 
spectral weights of the spinon and holon are redis- 
tributed with changing temperature, though the spinon 
and holon branches are not resolved in this experiment. 
The finite-temperature effect may also appear in the 
room-temperature ARPES data for SrCu02- 

Motivated by these ARPES measurements, we exam- 
ine the single-particle excitation spectra in the ID Hub- 
bard model at half-filling. We are particularly interested 
in the temperature dependence of the spinon and holon 
excitations at temperatures of the order of the magnetic 
exchange J w At 2 /U, where t is the hopping matrix el- 
ement and U is the onsite Coulomb repulsion. In the 
large-£7 limit, the single-partic le sp ectral weight was ob- 
tained in Refs. [a Q> IE S lii! IU LL-I- However, in the 
limit, the spinon excitations are not expected to exhibit 
temperature dependence. Therefore, such a simplified 
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picture of the large-J7 limit is not relevant for our discus- 
sion. 

In order to study the temperature dependence, we cal- 
culate the single-particle spectral weight of the ID half- 
filled Hubbard model by using the determinantal quan- 
tum Monte Carlo (QMC) and the maximum-entropy an- 
alytic continuation methods. For a complementary study 
to the QMC and a check of validity of our analytic contin- 
uation results, we also perform dynamical density matrix 
renormalization group (DDMRG) calculations for the 
spectral weig ht at zero temperature [Tl 111 El liH . 
It is noted that the single-particle spectrum of the ID 
Hubbard model has been previously studied by u sing th e 
QMC and the maximum-entropy techniques [Tj, [n| |2(| • 
In Ref. ^3 j the single-particle spectrum and density of 
states were calculated for U = At and T = 0.0625t at 
half- filling and at 1/6 doping. In addition, the velocities 
for the spin and charge excitations were obtained from 
the frequency and momentum dependences of the spin 
and charge susceptibilities. In Ref. [20|. the general fea- 
tures of the single-particle spectrum were discussed for 
U = 7.5* and T = 0.08t at half-filling. In these QMC 
data, however, the spinon and holon branches on the 
single-particle excitation spectrum were not resolved. 

The organization of this paper is as follows. In the 
next section, we show the U dependence of the spectral 
weight at zero temperature which was obtained by using 
the DDMRG method. In Section III, the temperature de- 
pendence of the spinon and holon excitations is discussed 
by combining the DDMRG and QMC results. The main 
purpose of this section is to resolve the spinon and holon 
excitations in the single-particle excitation spectrum of 
the ID Hubbard model at half-filling using the QMC and 
the maximum-entropy techniques. The discussion and 
summary are given in Section IV. 



II. ZERO-TEMPERATURE DYNAMICS 

The Hubbard Hamiltonian is given by 

H = -tJ](4 CT c l+ i, CT + H.c.) + [/J]n a n u , (1) 

i,er i 

where Ci^ a (cj CT ) annihilates (creates) an electron with 

spin a at lattice site i, rii = n^-f + n^i, n.i y(J — c| CT Ci j(T , t 
is the hopping integral along the chain axis, and U is the 
on-site Coulomb repulsion. We consider the half-filled 
case. 

In this section, we show the U depencence of the spec- 
tral weights of the spinon and holon branches at zero 
temperature. The results will be helpful for understand- 
ing the [/-dependence of the temperature evolution of 
the spinon and holon branches in the next section. The 
single-particle spectral weight is defined by 
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7T 



1 



H 



■ vy 



-Ck,1 



where Cfej is the momentum representation of the elec- 
tron operator with spin f, |0) and E are the ground 
state and the eigenenergy, respectively, and 7 is a small 
positive number. 

In the DDMRG method, the superblock is set to be 
64 sites, and we use the open boundary condition |l.l| . 
This is because the standard DMRG method gives the 
most precise numerical results when the open boundary 
condition is applied. However, there is a boundary ef- 
fect, which makes a photohole localized at the ends of 
the ID chain. For the single-particle spectral weight, it 
is not so easy to remove the boundary effect by increas- 
ing the system size. Thus, we add a potential of — trii 
to the ends in order to remove the boundary effect. It 
is noted that the spinon branch is almost unchanged by 
the potential, while the holon branch is shifted to higher 
binding-energy region by the order of O.lt. The momen- 
tum representation of the electron operator, Cfcj , in the 
open boundary system is defined by 



Cfe,T = 



— y^sin(fcZ)Q, 



(3) 
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in terms of an expansion in particle-in-a-box eigenstates, 
where k = nir/(L + 1) with n = 1,2,..., L, and L is 
the system size. In the infinite- chain limit, L — » 00, 
this transformation becomes equivalent to the standard 
Fourier transformation. It may be technically useful to 
note how to count the fermion sign when the state Ck,i |0) 
is generated, because the number of electrons at each site 
is neither 1 nor in the DDMRG bases of |0). Instead of 
cj f, we use an annihilation operator, c\^, which includes 
the fermion sign caused by the hopping from the Z-th site 
to one end of the system (environment) block |2l|. The 
density matrix we introduce is composed of four target 
states, i.e., the ground state ip a =i = |0), the final state 
after annihilation of an electron ip a =2 — Oct 1 0) , and two 
correction vectors ip a =3 = (Eq — uj a — H + £7) c/.,-)- |0) 
and -0q=4 = {Eo - u b - H + i-f)~ 1 Ck^ |0) being uu b = 
u> a +2j. The method is the so-called two correction vector 
method [15j . In the present work, the correction vectors 
are calculated by using the modified conjugate gradient 
method. By using these four target states, the density 
matrix p for each u> is defined by 
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where the indices i and j run over all bases of the sys- 
tem and environment blocks, respectively, ^2 a Pa = F 
and p a is set to be 0.25 for each a. We have checked 
that the results are almost unchanged for other sets of 
{p a }- The DDMRG bases are truncated from the eigen- 
states of this density matrix. Our method is equivalent to 
the variational principle for the mixed states |13|. Since 
p depends on u>, the truncated bases are optimized in- 
dependently for each u>. This is the reason why the 
renormalization works well for dynamical quantities. The 
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FIG. 1: (Color online) A(k, uj) versus ui for various values of 
U/t obtained by the DDMRG method. The momentum k is 
set to be 7r/65, which is the smallest value in the DDMRG 
calculations. The red line is the quasiparticle peak for U/t = 
0, obtained for k — by using periodic boundary conditions 
and a finite broadening of O.lt. The origin of frequency is 
shifted by ui s for the finite-(7 values, so that all of the peak 
positions occur at uj = — 2t. The value of the shift ui a is ~ 0.0, 
-1.134, and -2.80t for U = 4t, 8t, and 12t, respectively. 



truncation number which is sufficient for convergence de- 
pends also on uj. Here, the DDMRG bases are truncated 
up to to = 400. After convergence of the finite-system 
DMRG algorithm, we obtain A(k,uj a ) = — Im ( ip2 1^3 ) 
and A(k,uJb) = — Im (t/^IV^) /tt- At the same time, the 
spectral weights at an energy interval u> a < uj < u>b are 
calculated by using the same E and H. The broadening 
factor is set to be 7 = O.li, which is less than temper- 
atures treated in the next section. It is noted that the 
two correction vector method is efficient, but we need 
a higher truncation number than that in the standard 
DDMRG with a single correction vector. Then, we check 
the convergency for some uj points by using three target 
states with the single correction vector and by taking the 
truncation number up to m = 480. 

In Fig. ^ we show the U dependence of the spectral 
weight at k = 7r/65, which is the smallest momentum 
in the open boundary system. In finite-C/ cases, we find 
two peaks at the band edges. The peaks at low and high 
binding-energy sides correspond to the spinon and holon 
excitations, respectively, because both peak positions are 
equal to those predicted by the Bethe-anzats solutions. 
In this figure, these peaks do no show clear branch cuts 
due to the finite broadening factor of 7 = O.li. However, 
we have checked that the edge singularities at the peak 
positions become observable with smaller values of 7. 

In the case of U/t<8 (U/t>8), the peak height of 
the spinon (holon) branch is larger than that of the holon 
(spinon) branch. In the finite-t/ cases, the spectral shape 
of the holon branch is almost unchanged, while the weight 
of the spinon branch rapidly decreases with decreasing 



temperature. 



III. THE SPECTRAL WEIGHT AT FINITE 
TEMPERATURES 

Here, we discuss the temperature dependence of 
A{k, uj) of the ID Hubbard model. For this purpose, we 
present data obtained by using the determinantal QMC 
technique [2^ • With this method, we have calculated the 
single-particle Green's function 



G{k,T) 



J2^ kre (Tr^ +e Ar)cU0)), (5) 



where (• • •) denotes thermal averaging, Ci CT (r) = 
e Hr Ci a e~ HT , and T T is the Matsubara time-ordering op- 
erator. For temperature T (inverse temperature /3), the 
integral equation 



G(k,r) = 



+00 



duj ■ 



1 



(6) 



expresses G(k, t) in terms of the single-particle spectral 
weight 

A(k, uj) = Im G{k, iuj n — > uj + vy) (7) 

7T 



where 



G(k,iu) n ) 



/ dr e LUJnT G{k, t) 
Jo 



(8) 



and uj n = (2n+l)7rT is the fermion Matsubara frequency. 
We have obtained A(k, uj) from the QMC data on G(k, r) 
by solving Eq. © with the maximum-entropy analytic 
continuation method. 

We also present results on the static magnetic suscep- 
tibility at zero frequency, x(<?)i which is defined by 



X {q)= f drJ2e- iqe {m! +i (r)ml(0)), 
Jo 



(9) 



-,Ht 



where 777^ = — and mf (r) = e" ' m^e 

The determinantal QMC technique does not have the 
"fermion sign problem" |2^| at half-filling or in the ID 
case. However, the QMC algorithm which proceeds by 
using local updates of the Hubbard-Stratonovich spins 
produces results with long autocorrelation times at low 
temperatures and for large values of U/t j^. It has 
been noted that this is because the determinantal QMC 
algorithm with single-spin-flip moves does not explore 
the phase space in an ergodic manner for large values of 
U /t. In order to remove this problem, global Monte Carlo 
moves have been introduced, and this way the antiferro- 
magnetic structure factor for the 4x4 Hubbard lattice 
was calculated for U up to lQt |24|. However, the effects 
of global moves on dynamical quantities have not been 
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explored. Because of these reasons, we have restricted 
our QMC calculations, which use single-spin-flip moves, 
to parameter regimes where we did not observe ergodic- 
ity problems. In the following, we present QMC data on 
A(k,uj) for U = At at 0.25t < T < l.Ot and for U = 8t 
at 0.33t < T < l.Ot, and make comparisons with the 
DDMRG data obtained at T = 0. 

In obtaining A(k,Lo) from QMC data on G(k,r), we 
have used the maximum-entropy analytic continuation 
procedure described in Ref. [2|j| . As it is well known, the 
maximum entropy technique has finite resolution which 
decreases away from the Fermi level. Furthermore, at low 
temperatures and large values of U/t, the G(k,r) data 
exhibit long autocorrelation times [24|. In order to im- 
prove the accuracy of the maximum-entropy results for 
A(k,Lo), we have obtained QMC data on G(k,r) with 
good statistics. For example, the covariance matrix of 
G(k,r) used in the maximum-entropy technique always 
exhibited a continuous eigenvalue spectrum, as discussed 
by Ref. |2^|. In addition, the Bry an's and the classical 
maximum-entropy algorithms [25| produced similar re- 
sults for A(k,Lo). Furthermore, we have monitored the 
results of the maximum-entropy procedure for A(k,Lo) 
as the statistics of the QMC data on G(k, r) improved. 
These provide information about the reliability of the 
maximum-entropy images of A(k,Lo) presented in this 
section. 

In the following, we present QMC data for a 32-site 
chain with periodic boundary conditions, and make com- 
parisons with the DDMRG data obtained for a 64-site 
chain. Here, A(k,oo) is plotted in units of t — 1 , and 
A(k, lo) = A(ir — k, —u>) at half-filling. It is noted that the 
momenta which we can access in the DDMRG method 
are different from those in the QMC method, because 
different boundary conditions are taken into account 
in these methods in order to keep numerical precision. 
Then, the QMC result with momentum k is compared 
with the DDMRG one with momentum Lk/(L + 1), and 
the QMC result with k — is also compared with the 
DDMRG one with the smallest momentum tt/(L + 1). In 
order to make the difference small, we take the system 
size L as large as possible in the DDMRG method. 

Figure H shows the QMC and DDMRG results on 
A(k,uj) for U = At at half-filling. Here, the QMC results 
were obtained at T — 0.25i, and the DDMRG data are at 
T = 0. For k = 0, the DDMRG data show that the holon 
and spinon branches of A(k, lo) are at lo w — 3.5t and 
w — 2t, respectively. In addition, the insulating gap and 
small amount of spectral weight at lo ss 3.5t are observed. 
For k = 7r/2, the peak in the DDMRG results for A(k, lo) 
occurs at w — 0.7t. Hence, we deduce that A 0.7t for 
the Mott-Hubbard gap. As we go from k = to tt/8, the 
DDMRG results show that the spinon peak exhibits weak 
dispersion, while the holon peak moves rapidly towards 
the spinon branch. As k — > 7r/2, the spinon and the holon 
branches merge together. At T = 0.25t and for k — 0, 
the maximum-entropy images of A(k, lo) exhibit a peak 
at lo w — 2t and a shoulder centered at w —3.2t, which 
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FIG. 2: (Color online) Single-particle spectral weight A(k,u) 
for the ID Hubbard model at half-filling for U = 4t. The blue 
curves denote the QMC data at T — 0.25t for wavevectors 
< k < tt/2 for a 32-site chain with periodic boundary con- 
ditions. The red curves denote the DDMRG data at T = 
for wavevectors k = 0, it/8, it/4:, 3n/8 and it/2 for a 64-site 
chain with open boundary conditions. For definition of k in 
the DDMRG data, see the text. 



we attribute to the spinon and holon excitations, respec- 
tively. However, it is not possible to resolve the holon 
and spinon branches for k = tt/8 at T = 0.25i. We note 
that, at T = 0.25t, as k goes from 37r/8 to tt/2, the height 
of the peak in A(k, lo) decreases. This behavior is also 
observed for the DDMRG data at T = Q. For k = tt/2, 
the peak in A(k,Lo) at T = 0.25t is significantly rounded 
compared to the DDMRG data at T — 0, which includes 
an artificial broadening of 7 = O.lt. We attribute the 
rounding of the peak in the QMC data for A(k — tt/2, lo) 
to finite-temperature effects. For k = 3ir/8, we also ob- 
serve that, at T = 0.25t, there are single-particle exci- 
tations at lo ~ of which intensity decreases as k goes 
towards the zone center. 

Figure |2| shows the evolution of the maximum-entropy 
images of A(k,Lo) as the statistics of G(k,r) improve for 
U — At and T = 0.25i. This is useful in order to have 
an estimate of the resolution of the maximum-entropy 
method. We note that the dependence of the maximum- 
entropy results for A(k,Lo) on the statistics of the QMC 
data was also investigated in Ref. f° r the 12-site Hub- 
bard chain with U — At and T = 0.05t. In Fig. 01 A(k, lo) 
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FIG. 3: (Color online) Evolution of the maximum-entropy 
images for A(k, lu) as the statistics of the QMC data on 
G(k,r) improves. Here, A(k,oj) is plotted as the number of 
runs used for obtaining G(k, r) increases from 350 to 5000, 
where one run corresponds to one independent QMC simu- 
lation with 4 x 10 5 Monte Carlo sweeps. These results are 
shown for wavevectors k — it/2, it/4: and 0. 



FIG. 4: (Color online) Temperature dependence of A(k,u) 
versus lu for U = 4t and half- filling for wavevectors k = tt/2, 
7r/4 and 0. Here, the finite-temperature results are from 
maximum-entropy analytic continuation of QMC data for 
a 32-site chain with periodic boundary conditions, and the 
T = results were obtained by DDMRG for a 64-site chain 
with open boundary conditions. 



versus lu is plotted for k — tt/2, 7r/4 and 0, as the num- 
ber of runs increases from 350 to 5000. Here, each run 
corresponds to one independent QMC simulation with 
4 x 10 5 updates of the Hubbard-Stratonovich fields. We 
find that, for k — tt/2, the maximum-entropy results on 
A(k, lu) converge rapidly. However, for k = 0, the resolu- 
tion worsens, in particular, for structures in A(k, lu) away 
from the Fermi level. Hence, it is necessary to have QMC 
data with good statistics in order to resolve reliably the 
features away from the Fermi level. Here, we observe that 
the maximum-entropy images of A(k — 0, lu) for 2100 and 
5000 runs are similar. The QMC results on A(k, u>) were 
obtained by performing approximately 5000 runs. 

In Fig. 21 we show results on the T dependence of 
A(k,u) for U = At at half-filling for wavevectors k = 0, 
tt/4 and tt/2. At T = l.Ot > A, A(k = tt/2, lu) exhibits 
a broad peak centered at lu — 0, while at T = 0.5t A 
the insulating gap starts to develop. For k = 0, the T 
dependence is more involved. Here, at T = l.Ot £ A, 
a significant amount of the spectral weight is observed 
atwK — 3.5i, which corresponds to the location of the 
holon peak at T = 0. Furthermore, an additional peak 
at lu w and also spectral weight at lu > are observed. 



At T — l.Ot ^ A, there is a pseudogap at lu rj — 2t, 
which is the location of the quasiparticle peak for the 
nonintcracting system. As T decreases from l.Oi to 0, we 
observe that spectral weight from lu ks — 3.5< and from 
lu w are transferred to ui w — It to form the spinon 
peak. We think that the pseudogap at lu rj —2t is due to 
the onsite Coulomb repulsion and, for T > U, the peak 
in A(k = 0,lu) occurs at u> = —2t. However, we also 
observe that, for k = 0, there is significant amount of 
spectral weight at the frequency of the holon excitations 
already for T < U, while the spinon branch develops for 
T < l.Ot RJ J. The temperature evolution for k = tt/4 is 
similar to that for k — 0. However, the spinon and holon 
peaks are located closer, and it is not possible to resolve 
them at T = 0.25t. 

In Fig. we show results on the T dependence of 
A(k,u>) for U = 8t at half-filling. For k = tt/2, we ob- 
serve that the Mott-Hubbard gap is approximately 2At. 
In this case, we have the magnetic exchange J rj 0.5t, and 
we observe that A(k, lu) exhibits strong T-dependence. 
In particular, at T — 0.33< ^5 J and for k — 0, the 
maximum-entropy image of A(k, lu) shows a double peak 
structure of almost the same peak heights. The peak 
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FIG. 5: (Color online) Temperature dependence of A(k,u>) 
versus u for U — St and half-filling, plotted in the same way 
as in Fig. 4. 
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positions are nearly equal to those at T = 0, and thus 
the peaks at T = 0.33* are attributed to the spinon and 
holon branches. The spinon branch becomes observable 
at T = 0.33*. 

Finally, we show the magnetic susceptibility data. Fig- 
ures Efa) and (b) show QMC results on the temperature 
dependence of the magnetic susceptibility x(l) f° r U = 4t 
and 8*, respectively. In these figures, x{l) versus q is plot- 
ted for the same temperatures as those used in Figs. 0] 
and[S] For these U values, x(<? — * 7r ) increases rapidly as 
T decreases. Figure [7| shows QMC results on the temper- 
ature dependence of the uniform susceptibility x(q ~ * 0) 
for U = At and 8*. By determining the maximum of 
x{q — > 0), we estimate the temperature where the anti- 
ferromagnetic correlations develop. The temperature is 
estimated to be T w 0.5* for U = 4*, and T w 0.33* for 
17 = 8*. For [7 = 8*, the value 0.33* is close to J w 0.5*, 
while, for U = 4*, the value 0.5* is far below J w 1.0*. 



IV. SUMMARY AND DISCUSSION 



FIG. 6: Magnetic susceptibility of the ID Hubbard model at 
zero frequency x(l) versus q for (a) U — 4* and (b) (7 = 8*. 
Here, results are shown at half-filling for various tempera- 
tures. 
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Here, we summarize the numerical results obtained in 
the previous two sections, and discuss recent experiments 
in the light of them. In Section II, the U dependence of 
A(k ~ 0, oj) was calculated by the DDMRG method at 
T = 0. For the values of U presented in Fig. 1 (4* < U < 



FIG. 7: Temperature dependence of the uniform magnetic 
susceptibility x(l ~* 0) for U = 4* and 8t at half-filling. 
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124), the spinon and holon branches are clearly resolved 
for k 0. We have also checked that the spin-charge 
separation is observable for U down to l.Ot, though we 
do not show these data. In Fig. ^ the line shape of the 
holon branch is almost unchanged for 44 < U < 124, while 
the peak height of the spinon branch rapidly increases as 
U decreases from 124 to 44. 

In Section III, the T dependence of A(k,uj) was pre- 
sented at half-filling in the QMC method. A(k — > 0,w) 
clearly shows the double peak structure at T — 0.254 and 
0.334 for U — At and 84, respectively. At these temper- 
atures, the uniform magnetic susceptibility has a maxi- 
mum as a function of T, and the short-range antiferro- 
magnetic correlations develop. 

Let us first compare the numerical results with the 
ARPES data for Nao.9 6 V 2 05 5]. The value of U of this 
compound is calculated to be U ~ 12t in the finite- 
temperature exact-diagonalization calculations for the 
ID 4- J model 0- Therefore, the results for U = 84 may 
be applicable for comparisions at a qualitative level. The 
DDMRG (QMC) results are obtained at T = (0.334), 
whereas the temperature region where the ARPES mea- 
surements were carried out is estimated to be between 
and 0.334. In the numerical results on A{k, u>) for k = 0, 
7r/4 and 7r/2 shown in Fig. [S] the peak height of the 
spinon branch increases as T decreases from 0.334 to 0. 
In this temperature region, the threshold at low energy 
side is almost unchanged. For k — 0, the width of the 
lower Hubbard band becomes narrower with decreasing 
T, while for k = n/2, the width is almost independent 
of T for T < 0.54. These results are consistent with the 
ARPES data on Na .96V 2 O5. 

Next, we compare the results for U = 84 with the 
ARPES data for SrCu02 Q , since the appropriate value 



of U for this compound is estimated to be U w 104. 
In the ARPES data obtained at room temperature, the 
spinon and holon branches are resolved near k = 0. As k 
goes from to 7r/2, the spinon and holon branches merge 
into a single peak. The band widths of the spinon and 
holon branches are consistent with the numerical results 
shown in Fig. 5. The peak height of the holon branch is 
smaller than that of the spinon. On the pther hand, in 
the numerical results for A(k — > 0,oj) at T = and 0.334 
shown in Fig. the spinon and holon branches have al- 
most equal spectral weights. The discrepancy between 
the theory and the experiment may be attributed to the 
electron-phonon interaction. Within this context, it has 
been shown that the electron-phonon interaction broad- 
ens the holon branch more than the spinon branch |2fi| . 
The detailed origin of this remains to be clarified. 
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